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ABSTRACT 

Modal Analysis is becoming an integral part in the 
design phase and trouble-shooting aspects of noise and 
vibration related problems. Low cost digital hardware and 
increaSingly efficient software have facilitated the 
widespread use of this technique. 

A study iS made into the implementation of the impulse 
loading method for the purpose of determining the modal 
parameters of a structure uSing transfer function analysis. 

From this study, a practical modal analysis approach is 
outlined that is designed to be used without the aid of 
Specialized modal analysis units. Details of this procedure 
include: aspects of data aquisition, windowing in the time 
and frequency.domains, curve fitting and the use of curve 
fit values to estimate the modal parameters in question. 

Tests were performed on a continuous free-free beam to 
verify the technique by comparing the results with those 
predicted from theory. The technique was then applied to a 
practical engineering problem: dynamic analysis of an air 


compressor and its inertial base. 
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Nomenclature 


a, component of residue matrix 

vad Cara coefficients of numerator function, N(w) 

i geal residue matrix 

et Biol hei s coefficients of denominator function, D(w) 

[B(s)] system matrix 

S lower residual term of analytical curve fit 
function 

[c] damping matrix 

D upper reSidual term of analytical curve fit 
PUneT LON 

D(w) denominator function of Levy expansion 

i frequency 

q(t) time domain signal 

G(F) Fourier transform of g(t) 

[H(s)] transfer matrix 

eal identity matrix 

k pertaining to the kt® mode 

[K] stiffness matrix 

[M] mass matix 

n number of modes 

N number of time sample points 

N(w) numerator function of Levy expansion 

Px system pole 

Ss complex variable 

At sampling interval 

T length of time record 

{u,} modal vector 
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w(t) 


windowing function 


displacement vector 


complex constant 

amount of damping in exponential window 
coherence function 

error function 

damping ratio 

damping coefficient 


damped circular frequency 


denotes complex conjugate 
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1. INTRODUCTION 


1.1 F.E.M. and Modal Testing 

Exact, closed form solutions describing the dynamic 
behavior of a body exist for but a few simple homogeneous 
configurations. For more complex structures, it is necessary 
to use alternate methods to determine their vibrational 
characteristics. Commonly employed procedures are finite 
element analysis and modal testing methods. 

In finite element analysis, the structure is divided 
into a finite number of discrete elements. The equations of 
motion of the individual elements along with the 
constraining boundary conditions, together describe the 
overall motion of the body. Once this mathematical model has 
been formed, it can be used very effectively to study both 
the static and dynamic behavior of the body due to any 
number of different loading conditions. The effects of 
changes in design and material selection can at once be 
determined, whereas actual prototype testing would be a much 
Slower and more expensive route to choose. 

Designing using the finite element method, does have 
its problems however. They include; the large computer 
programs and memory space required to model the structure 
adequately and the often large discrepencies between 
predicted and actual test results. Direct dynamic testing, 
then is necessary to determine the validity of any computer 


modelling that is carried out. 
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Dynamic or modal testing is normally of one of two 
types; normal mode or transfer function methods. In normal 
mode testing, resonant modes are excited one at a time 
(using shaker tables) and the natural frequencies, damping 
and mode shapes can be determined. This process however, is 
difficult to use when two adjacent modes are close together 
because of the problem of exciting only one mode at a time, 
especially if the damping is appreciable. Another problem 
with this method is its limited versatility. Field testing 
1s very impractical because of the large amount of equipment 
to set up, and the analysis of large structures is limited 
by the size and number of shaker tables available. 

The transfer function method is superior to the normal 
mode method, in that all of the dominant modes can be 
measured simultaneously. As well a large variety of single 
unit exciting methods can be employed to force the structure 
which helps keep equipment and time requirements to a 
minimum. By using the Fourier transform to determine the 
transfer function between various points on the tae 
structure, and curve fitting these experimental transfer 
functions to a standard form, it is possible to quickly 


determine the vibrational properties of the structure. 


1.2 Practical Modal Analysis 
Transfer function test signals are usually digitized 
rrgnt on locatron using duel’ channe? FFT) analyzers For the 


purpose of this study, however, it was deemed more useful to 
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record raw field data on a portable multichannel, FM tape 
recorder. The main reason for choosing this route was that 
many applications of the type of analysis require the taking 
of data in dirty or hazardous environments (e.g. chemical 
refineries). In this way, test data could be quickly stored 
on tape and later played back for analysis. A second and 
perhaps aS important a reason as the first, is that it 
allows flexibility in the processing of the signals. Tests 
can be done on various frequency ranges, all using portions 
of the same test data. Unless special provisions are made, 
(i.e. using high sampling rates and large sampling times) 
purely digital test systems require additional tests if 
different frequency ranges are to be used. Another 
flexibility that this arrangment lends itself to is the fact 
that a multichannel FFT unit is not required to process the 
information. Instead, the analog data may be filtered and 
digitized with conventional equipment. 

Privately developed or commercially available software 
can then be employed to reduce the data to yield the desired 
modal parameters. This approach allows a company or 
consulting firm to be able to do complete modal analysis 
with equipment (FM recorder, filters, digitizers and 
computing facility) that is usually already available. The 
complete modal analysis units on the other hand, tie up a 
lot of capital in a system that in many cases has limited 


flexibility and a short half-life before becoming obselete. 
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1.3 Research Outline 

The intent of this study was to investigate the 
flexiblility of implementing a practical approach to modal 
analysis. This was to include both the writing of the data 
acquisition and modal analysis software as well as a 
description of what type of equipment is required to 
successfully carry out this procedure. 

Chapter two of this report contains the formulation of 
transfer functions in the Laplace domain from the equations 
of motion for a structure. It also shows how results 
obtained from a Fourier transform of a time record relate to 
the Laplace domain representation of the body. This chapter 
concludes with a discussion on the limitations and 
assumptions of the discrete Fourier transform. 

Chapter three deals with aspects of curve fitting the 
experimental data to the analytical representation of the 
transfer functions that were developed in the previous 
chapter. It shows how an approximate linear method can be 
used to generate good initial estimates for the modal 
parameters in question and how a simultaneous, multi-degree 
of freedom curve fit can be used to accurately complete the 
process. This is followed by a modal vector identification 
scheme that can deal with either normal or complex mode 
Shapes. 

The following section, Chapter four, addresses some of 
the practical considerations that one must be concerned with 


when using transfer function analysis. The first deals with 
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the types of excitation procedures commonly used and why 
transient testing was chosen in this study. This is followed 
by a discussion into the kind of equipment that is required 
to first of all record the field data, and later play back 
the signals to be digitized and analyzed. The final part of 
this chapter looks into the coherence function and its use 
in determining some properties about the quality of the 
generated signals. 

Chapter five discusses the set-up and results of two 
experimental tests. The first test was done on a free-free 
beam to allow easy comparison between theoretically 
predicted and experimentally derived results. The second set 
of tests were done on an air compressor base as an example 
of a practical engineering application of the procedure. 

The final chapter lists the conclusions to the study 
and points out possible areas of future research that stem 


from this project. 
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2. THEORETICAL CONSIDERATIONS 
To understand the transfer function method of testing, this 
chapter begins with the development of the transfer 
functions from the equations of motion. It is then shown how 
certain properties of the set of transfer functions for the 
Structure, allow greatly simplified testing techniques when 


employing this method. 


2.1 Transfer Functions in the Laplace Domain 

Consider a continuous structure whose dynamic behavior 
is to be determined. It will be modelled by a series of 
discrete masses interconnected with linear springs and 
viscous dampers. The motion of this configuration may then 
be described by a set of second order linear differential 
equations with constant coefficients. In the time domain, 


these equations would have the form: 


Mince) tart kG rocker. jak | KoleeCh) haat Gt 7 (27 hr) 
where: 1£(t) 3 s=Japplied foreing ‘vector (nx 1) 
{x(t)} - resulting displacement vector (nx1) 
{x(t)} & {¥(t)} - first and second derivatives of 


{x(t)} with respect to time (nx1) 


[Mii Gis bk] mass, damping and stiffness 


matrices (nxn) 


Although it is not necessary to do so, we will only 


consider the case where [M],[C] and [K] are symmetric 
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matrices with real valued elements (see Potter, [13]'). 
Taking the Laplace transform of equation (2.1) and 
considering the case where the initial velocities and 


displacements are zero, at t=0, it will follow that: 
[M]{x(s)}s*.+{C]{x(s)}s +[K]{x(s)} = {£(s)} G22) 
Defining the system matrix as, 
EBCS)ifiee [Mls cm + LC ]is tok) 253°) 
equation (2.2) can be simplified to: 
[B(s)lix(s)} = {f(s)} (2.4) 
The transfer matrix of the system, [H(s)], is defined as: 
[HCs)Ja= (BCs) ys ¢ 2), 


Considering equation (2.3) and noting that, 


[H(s)] = adj[B(s) ] 
det [B(s)] 


it is apparent that each element of [H(s)] is a rational 
function of the Laplace variable, 's'. Therefore each term 
in the transfer matrix can be expressed in partial fraction 
form. The zeroes of the det[B(s)] (the poles of the transfer 
matrix) always occur in complex conjugate pairs except when 
the system is critically damped or over-damped (poles will 


be real valued in these situations). 
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Assuming the poles are of 


be written as, 


unit multiplicity, [H(s)] may 


S-p, S-pi* S-Dn S-Dn* 
Or more simply as, 
n 
(H(s)] = = j[a,J + [A,*] (2561) 
k=1/S-p, S-Dy* 
where: [A,] - complex residue matrix, (nxn) 
Sea , th pole of the system 
n - total number of modes 
"*' = denotes complex conjugate 


2.2 Residue Matrix Properties 
To determine some useful properties of the residue 
matrix [A,], equation (2.6) is pre-multiplied by the system 


matrix [B(s)] and the scaler (s-p,). 


(s-p,)[B(s) ][H(s)] = aera 5 ss (OAR eet cu 


S-p.* 


From equation (2.5) we see that [B(s)][H(s)] = [I] (the 


identity matrix), therefore the left hand side of the above 
equation reduces to (s-p,)[I]. If we let s=p,, the left hand 
side will be identically zero and the right hand side will 


be reduced to [B(p,)][A,], therefore: 
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[B(p.)][A.] = 0 (20) 


In a Similar manner, post-multiplying equation (2.6) by 


[B(s)](s-p,) and allowing s=p,, it will be seen that 
[A, ]J([B(p,)] = 0 G25.6)) 


Now consider equation (2.4) with the forcing function, 


{f£(s)}, set equal to zero. 
[Sts jtseCsef =.10 


These homogeneous equations correspond to the free, damped 
vibration of the system. Non-trivial solutions will only 
exist if the det[B(s)] = 0. For a value of 's' that 
satisfies det[B(s)] = 0, say s=p,, there will also exist a 


unique solution vector, {u,} such that: 
BE Cow fu 630 (2.9) 


{u,} is known as the modal vector or mode shape for the pole 
S=p,. 

If we now consider equations. (2.7) & (2.8) in the light 
of equation (2.9), it can be seen that both the columns and 
the rows of the residue matrix [A,], must be linear 
combinations of the modal vector {u,}. [A,] can then be 
replaced by, 

fA) sSab ltt, bias, bPep esta, fa, Pes (2810) 
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10 
where 'a,' 1S a complex constant. In general the elements of 
the modal vector {u,} are complex. If however, the system is 
proportionally damped?, the modal vectors will be purely 


imaginary. Substituting equation (2.10) into equation (2.6) 


yields, 


[H(s)] = ; Ot 33 er fueelins)) (ee) 
k=] SD; Seis 

Due to the properties of the residue matrix (equation 
(2.10)), it is only necessary to make measurements that will 
complete one row or column of the matrix to completely 
define the system. This can be accomplished by exciting the 
Structure at one point and recording the response at all 
other points of interest (simultaneously or one at a time). 
Alternatively the response can be measured at one fixed 
point and then allow the forcing ada to be positioned 


at various points on the structure. 


2.3 Frequency Response Function 
Between any two points on the test structure, the 
transfer function can be found by forming the following 


operation, 


where X(s),; is the Laplace transform of the response at a 


point 'i' due to the forcing function and its corresponding 


2For a more general discussion, see Caughey Sar 
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For the purpose of this study, the complex variable 


"Ss, car be written ‘as, 
S=oat Jw 


’ 


where 'o' is the damping coefficient and 'w' is the damped 


angular frequency. A complex pole of a system is written as 
Dio = "Ores JO. 


The damping coefficient of a pole, o,, iS negative for a 
Stable system, and the natural frequency, o,, 1S actually 
the damped natural frequency for the particular mode. The 
undamped natural frequency (resonant frequency) is found 


from 
Biya es Cae tu ips in ee 
and the damping ratio is determined from 


ba 5s ON 
Wkn 
Both the real and imaginary components of the transfer 

function, H(s), map out continuous surfaces over the Laplace 
plane (except at pole locations). To be able to fully 
describe the surfaces and represent them in the form of 
equation (2.11) it is mot necessary to make measurements 
over the entire s-plane, because the surface is described by 


an analytic function. 
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In practice, the transfer functions are measured along 
the imaginary axis. What is in fact done is to evaluate a 
Fourier transform rather than a complete Laplace 
transformation. This is equivalent to letting s=jw (i.e. 
o=0) and computing the transfer function along this line. 
This representation of the transfer function is known as the 
"frequency response function'. From these measured values, 


the entire s-plane representation may be formed if desired. 


2.4 Discrete Fourier Transform 


The Fourier transform is defined as; 


G(F)) = J -alt)exp(qj2xft) dt. 
—o 
Practical limitations necessitate that the signal record be 


of finite duration, say from t=0 to t=T, such that; 
T 
Git i= lege) exp] j278tar:. 
0 


The analog signals, that are to be analyzed, are normally 
digitized to allow efficient handling and processing of the 
data. This discretization, coupled with the time limited 
Sample record, reduces the integral transform to a discrete 
Fourier transform of the form; 

N=] 


GCA eat Nee Deg ita exp(-j2atee, } (0041.2)) 
n=0 
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where 'N' is the number of samples of the time domain 
Signal. 

The 'windowing' in the time domain has the effect of 
limiting the maximum obtainable frequency resolution (Af = 
1/T). Any increase in frequency resolution for a given test 
must be accompanied with longer sampling times. The finite 
number of time sample points recorded (typically N=1024), 
limits the maximum obtainable frequency, fm, that may be 


determined. The sampling period At, is given by 


ancery, 15 found Erom 
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2A aes 


To obtain information in higher frequency ranges, 
shorter sampling periods must be employed. If N is to stay 
fixed however, this dictates lower (and often unacceptable) 
frequency pac cuaaonee To alleviate this problem, 'zoom' 
transformations are used, which allow the operator to choose 
the upper and lower frequencies over which the Fourier 
transformation is to be performed. This then allows the 
available number of spectral lines to be used in a specific 
bandwidth of interest rather than from DC to some maximum 


frequency (see [12]). 
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2.5 Windowing 

The digitizing of the time domain signals are subject 
to Shannon's sampling theorem which states that in order to 
avoid aliasing or 'folding back' in the frequency domain, 
the sampling rate must be greater than twice the highest 
frequency content of the waveform to be analyzed. By passing 
the signals through low-pass filters however, does not 
ensure that the frequency content above some designated 
frequency has been attenuated to zero. It would be necessary 
to sample at an infinite rate to avoid any aliasing for a 
general waveform. Sampling at a finite rate then, will 
involve folding back to some degree (except for a band 
limited function). Its effect however, can be safely ignored 
if the frequency content of the signal that is folded back 
is small compared to that in the frequency range of 
interest. 

In typical waveform analysis, the analog signals are 
first of all passed through sharp roll-off, low-pass filters 
(the cut-off is set at the highest frequency of interest, 
fm). These waveforms are then digitized at a rate of f,, 
where f,>2fm. Any frequency content of the signal greater 
than f,-fm will be folded back to the frequency range in 
question. However if f,=2.56f, and the low pass filter has a 
roll-off rate of 96dB/octave, then the signal folded back 
onto the frequency f,, would be attenuated by 62dB. This 
degree of suppression is normally more than adequate to take 


care of any potential aliasing problems. 
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Once the signals have been properly filtered and 
digitized, a discrete Fourier transform is performed on 
them. Examining equation (2.12) it is seen that the 
frequency domain, in addition to the time domain signal, is 
considered in a discrete rather than in a continuous sense. 
This 1s necessary because for digital computation, only 
discrete values of the frequency response function can be 
formed. It is shown in [3], that this ‘uniform sampling' in 
the frequency domain implies that the time domain signal be 
periodic. This means that any time domain signal that is to 
be analyzed, must be representative of one period (or an 
integer multiple of it). If this is not true, significant 
distortion of the frequency spectrum may result. To reduce 
this effect, time domain windows are used, which force the 
Signal to be periodic within the sample time. This is done 
by simply multiplying the measured, digitized signal by a 
function that tapers off to zero at either end of the 
window. A commonly used window in the study of continuous 


Signals is the Hanning window, which has the functional form 
wt) s=-05500 —scos(27t/T)} 0Stst. 


Although this, or any other applied window, does not totally 
eliminate the distortion (the windowing process itself adds 
some predictable amount of distortion to the signal) caused 
by a finite time sample, the windowing reduces the 
distortion to a level from which meaningful results can be 


derived. 
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During the course of this study, impulse testing was 
employed aS an exciting force in the determination of the 
transfer functions for a structure. This type of excitation 
resulted in transient waveforms for both the forcing and 
response signals. Typical experimental signals, that have 
been low-pass filtered, are shown in figures 2.1 and 2.2. 
When computing the Fourier spectrum, these transient 
waveforms require different windowing functions than those 
used with continuous Signals. In fact, the input force 
Signal does not even require a window because it is 
self-windowing (that is, it occurs completely within the 
allotted time block). The response signal on the other hand, 
which is composed of decaying sinusoids, may not have died 
out completely in the time record and must therefore be 
windowed. It is apparent that applying a Hanning window to 
this type of signal would not be appropriate, because in the 
process of forcing the signal to zero at both ends of the 
record, 1ti as, 1.0St,; +the, majority: of ats, information that 
was Situated at the beginning of the time record. For this 
type of signal, an exponential window is used. A commonly 
used form of this window weights the data by a factor of one 
at the start of the window to a factor of 0.05 at the end of 


the time block. This can be expressed as, 
wttyi= exp(-Bt} 


where B = 10g(0.05)/T. It should be noted that unlike the 


Hanning window, the exponential window does not force the 
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Figure 2.1 Forcing Signal 


Figure 2.2 Response Signal 
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Signal to be zero at both ends of the time record. The 
left-hand edge of the response signal is already zero due to 
the fact that the structure has not yet been excited, and so 
there is no need for weighting at this point. As the 
response signal continues in time, it decays exponentially 
and the effect of the window is to increase the apparent 
damping in the system, by bringing the right-hand edge of 
the signal as close to zero as desired. When the modal 
parameters are determined from the exponentially windowed 
response functions, the measured damping values (0) will be 
the sum of the actual damping present in a particular mode 
and the amount of damping present in the window itself (f). 
To obtain the true damping for a mode then, it is only 
necesSary to subtract the damping induced by the window from 
the total measured damping value. It should be noted that 
the apparent increase in damping induced by the exponential 
window, can cause tightly spaced, heavily damped modes to 
become very difficult to analyze and therefore may require 
additional tests with increased frequency resolution. 

The exponential window has an additional feature which 
makes it even more useful as a windowing function of the 
response signals. This is due to the fact, that as the 
response signal becomes smaller at the end of the time 
record, the signal to noise ratio is significantly reduced 
(especially in heavily damped systems). This would tend to 
add unwanted noise into the signal. However, when the 


exponential window is applied, it most heavily weights the 
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Signal at the beginning of the time block where the signal 
to noise ratio is highest, while the potentially noisy 
portion of the signal is largely suppressed. 

As stated earlier, no windowing is necessary for the 
impulse forcing function to meet the ‘periodic time domain' 
criteria of the discrete Fourier transform. Windowing is 
used on this type of signal however, to reduce the inherent 
noise in the Signal that is due to the short duration of the 
impulse compared to the total sample length. A simple type 
of eda that can be used is a rectangular function. With 
this window, the actual impulse portion of the signal, is 
given a weighting of unity and the data beyond that point is 
set to zero so that any stray noise in that section of the 


Signal 1S removed. 
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3. IDENTIFICATION OF MODAL PARAMETERS 

Once One “row Or CoOLlummmor the transfer “function matrix 
has been measured, it is then necessary to determine the 
modal parameters (natural frequency, damping coefficient and 
mode shape) for each mode in the frequency range of 
interest. Assuming a linear model is valid for the test 
Structure (as in equation (2.1)), the corresponding transfer 
functions may be represented by a simple sum of the response 
curves for each of the natural frequencies of the system. 

Figure 3.1 shows an example of a transfer function 
(solid line) and its contributing single degree of freedom 
components. Notice that the net transfer function lies below 
some of its contributing components. At first this would 
appear to be an inaccurate summation. It must be remembered 
however, that the transfer function is complex valued and 
that the overall PQeosioe Bupeeies is the sum of the real 
and imaginary parts of the single degree of freedom systems 
(only the magnitude is shown in figure 3.1). 

- The density of the modes will dictate what type of 
analysis is necessary for extraction of these parameters. 
When the natural frequencies of the system are well spread 
apart and the damping is light, single mode methods are used 
which make use of the data in the immediate area of the 
resonance peak. On the other hand, a multiple-mode technique 
is required to separate and identify a system with tightly 
spaced peaks and large damping factors. Here it is necessary 


that all of the modal parameters be identified 
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Figure 3.1 MDOF Transfer Function 


simultaneously. The methods employed are often iterative in 
nature and the initial values for the parameters to be used 
in the computations are typically obtained from single mode 
approximations. | 
In modal analysis, the process of identification ofea 

system is often termed "curve fitting'. This is because 
experimental data is frequently fit (in a least-squares 
sense) to an analytical transfer function. A common function 


used to fit the data is shown below. 


nN 
H( jw) = 53 a a (3.1) 
=1| 7 
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This equation is obtained by substituting s=jw into equation 
(2.6) (as the transformed experimental data is in the 
frequency domain) and by considering one element in the 
residue matrix [A,] (and the corresponding element in 
[A,*]). The non-linear nature of this function often makes 
it difficult to successfully carry out the identification 


schemes. 


3.1 Levy Curve Fit 

Generally, single mode methods are the simplest to 
implement due to approximations that can be made in 
estimating the parameters of equation (3.1). These 
techniques however are slow and often do not make use of all 
of the available data around the mode of interest to obtain 
a best possible estimate of the parameters. In addition to 
this, it was desired to employ a method that would work 
equally well in situations that would warrant either a 
Single or multiple mode analysis. 

To satisfy the above criteria, a complex curve fitting 
scheme devised by E.C. Levy was chosen to be used. (A short 
discussion outlining the major steps in the development of 
this technique will be presented, however an interested 
reader is encouraged to consult the original paper outlining 
the method [10].) 

Levy's method rewrites the transfer function of 
equation (3.1) as a ratio of two frequency dependent complex 


polynomials. This takes the form, 


~. 7h? : 7 —_— 7 
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Or more Simply as, 


G(jw) = N(w). G33) 


Let the actual experimental data at frequency ‘w,', be given 
as H(jw,). The error at this point is defined as: 

EXO) "H.C jo.) — CUjOp) eb jay)” =" Na) (34) 

D(w, ) 

Utilizing conventional least-squares analysis on the error, 
e(w,), this would yield a system of simultaneous, nonlinear 
equations to be solved. Iterative techniques are generally 
involved which require good starting values and some degree 
of operator experience. To avoid this problem, the error is 
multiplied through by the denominator of the analytical 


function. This results in: 
Ey )F MRE Dee) =H Os) Dope i LO ale, (3n5)) 


which is now a weighted error function. When the method of 
least-squares is applied to the modulus squared of this 
function, it results in a set of linear algebraic equations 
that.can, be e€asay solved for the coeffiictents Ac. 4. Aw, By 
vu Bo. (Of, equation .(3.2)). While. thisutype, ot approach. at 
first glance looks very promising, it produces 
unsatisfactory answers when large frequency ranges are 


encountered or when the denominator term, D(w), of the 
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analytical function has a wide range of values (like that 
which occurs near a pole). This is due to the weighting of 
the error function by the denominator term. It has been 
shown by Sanathanan and Koerner[20] that these effects can 
be virtually eliminated by a simple iteration to eliminate 
the weighting. The error function of equation (3.5) is 


modified to the following form: 


e'"(@,) # €' (oy) = elu, )D(O, ); G34/6)) 
D(w,) i -1 D(w,)i-s 


where the subscript '1' corresponds to the iteration number. 
Because D(w,) is not initially known, it is set equal to one 
for the first iteration (i.e. the first iteration is 
identically the weighted least-squares minimization). 
Ensuing iterations converge very rapidly. Five to ten 
iterations are normally required to have the parameters 
agreeing to within 0.1% from one iteration to the next. 

This technique proved to be very successful due to its 
speed, ease of application and agreement with the 
experimental data that it was meant to fit. The extraction 
of modal parameters with tightly spaced, adjacent modes were 
also correctly identified using this simultaneous curve fit. 
However, it strayed from the true parameter values when it 
was applied to fit two or more peaks simultaneously that 
were separated by a large number of data points. This was 
not found to be a problem because when the peaks were widely 
separated, they could be safely treated as single degree of 


freedom systems with little interaction from other modes. A 
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Simultaneous fit was therefore not necessary. 

When applying Levy's method, it must be decided how 
many terms of the numerator and denominator will be used in 
any particular curve fit. (see equation (3.2)) At first, one 
might select a large number of 'A' and 'B' coefficients and 
let the ensuing calculations determine how many of the terms 
were necessary (i.e. the terms that were not used in the fit 
would equal zero (or nearly so)). This 'blind' approach can 
lead to problems, because when many excessive terms are 
used, ill-conditioned numerical problems are frequently 
encountered. This problem need not occur however, because it 
is very straight forward to calculate the exact number of 
terms that should be used in the fit. For a single degree of 
freedom, the transfer function takes the form: 


H(jw) = _a tev atnen C357) 
Jo-p jJw-p* 


where: a ave resjay 


To + Jwo 


O 
i) 


Substituting the expressions for 'a'and 'p' and putting 
everything over a common denominator, equation (3.7) can be 
written as: 


H(joym= =2aro> t+ aywo) + (2a) jw (3)6)) 


os -Btta5* - (200) for] o- 


Dividing both numerator and denominator by ‘007 + wo*", 


equation (3.8) could be written in the form 
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H( jw) = Ao + A, (jw) 
WFP Shia) + Be ( jw)? 

where Ajo, A;,, B; and Bz are real constants that would be 
determined from the Levy curve fit. For each additional mode 
that is to be fit, it can be shown that it is necessary to 
add two additional terms to both the numerator and 
denominator. So that for two simultaneous modes, for 
example, it 1S necessary to determine the constants Aj, to A; 
ands B,-tovB,. 

Once the Levy coefficients have been determined, it is 
necessary to find the poles of the system, (i.e. to convert 
from the: form of equation (3.2) “back to that of “equation 
(3.1). This is done by calculating the complex zeroes of the 
analytical denominator function (e.g. the zeroes of '1 + 
Baja) i+ Be (ju)? + B,(ja)*+ 8, (jo)* " for the two degree of 
freedom case). These zeroes typically ks in complex 
conjugate pairs, but they can be purely imaginary if the 
damping is zero. It is then a simple matter of performing a 
partial fraction expansion (using the computed poles) to 
finish the transformation. (Appendix A lists a Fortran 
program that computes a Levy curve fit and then converts the 


derived parameters to the partial fraction form.) 


3.2 A True MDOF Curve Fit 
If the curve fitting process hasPnot includedwall of 
the modes of the structure (always the case for a continuous 


system) the curve fit will be inaccurate to some degree. 
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This is because the "talls' of the modes outside the 
frequency range of interest, will have been measured during 
the data acquisition process, (there is no way to filter 
them out) but no terms were included in the curve fit 
procedure to account for them. To choose what type of 
function would best suit the tails of the modes outside the 
desired range, it 1S necessary to examine the transfer 
function for one mode. Rewriting equation (3.9) again: 


H(jw) = Ao + A;( jw) 
1s B, (jw) au B2( jw) * 


it is seen for small values of 'w', that the transfer 
function may be represented by a real constant (in this case 
"Ao'). For large ‘'w' values on the other hand, the transfer 
function is proportional to 1/w* (again real valued; note 
that the constant 'A,' is zero when proportional damping is 
present or nearly zero when light, non-proportional damping 
exists). 

If it is desired for example, to fit the modes ina 
certain frequency range (w, to w,) of a transfer function 
and there exists other modes outside this range, the 
following analytical function would best represent what is 


occurring inside the range of interest. 


Jem. 
HC jo) eGR ie ay He an ae (3.10) 
aw? k=k; | Jw-px JO-p. * 


where 'w' is restricted to w,Swsw, and the modes from ‘'k;' 
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would take care of the effects of the tails of any modes 
below 'a,, and the constant term, ’D’', would similarly 
compensate for those modes above '‘'w,'. 

With this analytical model to work with, a method must 
be devised to correctly identify the constants in the 
expression. Ideally all of the parameters of equation (3.10) 
should be identified simultaneously and preferably using a 
least-Squares technique to smooth out the effects of any 
random error. Because of the nonlinear nature of this 
problem, good starting values (initial guesses) are 
necessary to increase the likelihood of convergence in the 
iteration scheme. As previously noted, Levy's method gave a 
good approximation to the modes within the range without 
including the effects of the modes outside the range (i.e. 
the term 'C/w* + D'), and therefore would give good starting 
values for the terms 'a,' and 'p,' (remembering that the 
subscript 'k' indicates the kt® mode in the range of 
interest)<« To initializes 'G'.andew'D'; thes difference-isi now 
formed between the experimental data and the Levy curve fit 
model. The real part of this complex difference is curve fit 
torpthes functions’ C/ose+aD*.) Althoughathise@fstamayynotegive 
accurate estimates of 'C' and 'D', it at least gives 
reasonable starting values to the two parameters. 

All of the parameters of equation (3.10) have now been 
assigned starting values for use in the simultaneous curve 
fit. The error function, e, that is to’ be minimized using 


the least-squares technique, is normally defined in the 
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following way: 


t 
PS ee ommnr HexD (jae) te 
m=s 
where: H(jwm) - analytical transfer function 
(e.g. @€quation (3. 10)) 
Hexp(jwn) - experimental transfer function 
Ol estesa | en tnd icacesmmodusus 


The minimization can now take one of several routes+7 1 
can proceed in a purely iterative fashion with all of the 
parameters optimized this way or, the minimization can go 
back and forth between solving for the parameters ‘'a,', 'C' 
and 'D' in a linear fashion and for the parameters 'p,' with 
an iterative approach. If some values of the parameters are 
known to be accurately determined, these can be held fixed 
while minimizing with respect to the other variables. 

In this study, the first of these three methods was 
chosen due to the ease in programming, that 1s, performing 
the minimization of all of the parameters with an iterative 
approach. This method of computation was fairly costly in 
terms of computer time and it is suggested that the second 
and third methods could be employed to speed up the process. 

The following three figures show the typical progress 
of a curve fit using the above outlined procedure. Figure 
3.2 shows the results of the initial Levy approximation (the 
solid line indicates the experimental data). The transfer 


function curves are only shown in the frequency range cf the 
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curve fit because it is only in this range that the curve 
fit is Walid: Figure 3.3 is the same as figure 3.2 except 
that it now includes the preliminary estimates of the 
coefficients in the term 'C/w? + D'. Notice the apparent 
"vast improvement' in the fit, just due to the addition of 
this term. At this point, the fit is not expected to be 
perfect, even though all of the terms of the analytical 
function (equation (3.10)) have been assigned values. This 
is because these coefficients are estimates only, and were 
not calculated simultaneously 

The final stage of the curve fitting procedure is the 
Simultaneous least-squares minimization of the error 
function, equation (3.11), with respect to the parameters 
Malet) Dy ip (Ges and! DD’... ASi can be seen in figure 3.4, the 
analytical curve fit falls very closely beside its 
experimental counterpart. 

The largest differences between the two occur where the 
curves reach their minimum values. At these locations, the 
magnitude of the error is small (though it appears large 
with the log scaling) and therefore has little effect on the 
Simultaneous curve fit. If an improved fit is desired, 
weighting would be necessary to emphasize the errors in 


these regions of the graph. 
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3.3 Mode Shape Determination 

The final curve fit parameters of equation (3.10) 
describe an experimental transfer function. The values 
obtained however, offer little insight into the problem at 
hand unless they are compared to, and used with, other 
transfer functions obtained from different testing points. 
The relationship between the transfer function parameters 
will become clearer in the following discussion. 

As noted in the previous chapter, the properties of the 
residue matrix [A,] made it necessary to only determine one 
row or column of the transfer matrix [H(s)]. Rewriting 


equation (2.4) as: 
fxd sdje eal Bs dei fics } 

and substituting equation (2.5), results in 
{x(s)} = [H(s)]{F(s)} 


To determine one row of the transfer matrix it would be 
necessary to measure the response at one point on the 
structure, due to forcing functions being applied, one at a 
time, to all of the other points of interest on the 


Srriictei ness tiie: 


H,, from response at point 1 due to force at point 1 
H,2 from response at point 1 due to force at point 2 


etc. 


In a Similar manner, one column of the transfer matrix would 
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be determined by forcing the structure at one point and 
measuring the response at all other desired points on the 
Structure (simultaneously or one at a time). 

In determining one row of the transfer matrix, for 
example, one row of each of the residue matrices have 
necessarily been determined. Consider a two degree of 


freedom system. 


Cue) Ci) (2) (2 *) Gt). Gtr) G25) (2 *) 
Ai1 tay; tay; tay Ai2 +ta;2 +Aa;2 +A;2 
HUS)G =) (S-D; S—p1* StDe S=ps*). Spi) S=Di*1S—-p2_S-p2* 
H2; H22 


where the ‘'a' and 'p' values are complex constants 


determined from the curvefit. This can be written as 


Gale) (aim) (1 *¥) (1 *) 
ai1 ai2 aii a12 
H(s) = ? ? + ? ? + 


S-Di S-pi* 


(35018) 


where the symbol '?' represents those values that were not 
determined from testing. Using equation (2.10), the first 


residue matrix can be represented by 


[A,] = Ja,u,u, @,U,U2 C3e12)) 
@,yU2U;, @1U2U2 


Matching common terms from equations (3.11) and (3.12), it 


is seen that 
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Cores 
Ari = @1U,U;, and 


CT 
C425 |S) Oy 2:, 
The components of the modal vector u;, uz are relative 
Quantities and may be calculated any number of ways (a; is a 
complex constant and u;, uz are the complex mode shape 


Gab) 
vectors). Dividing a,, by itself, yields a real number of 


unity value. Dividing same by atyety produces a complex 
mumber Sthat-ista *ratiosof*ao/7u, .oThe \form.ofxsthis. ratio 
Specifies the nature of relationship between uz and u,. If 
the ratio is real valued, a normal mode is present. The sign 
of this real number, positive or negative, indicates whether 
U2 1S in phase or out of phase with u,, respectively. If the 
ratio is zero, it indicates that the test point is a nodal 
point for the mode in question. If forming this ratio 
results in a complex number, it indicates that a complex 
mode is present. A complex mode indicates that when the 
Structure is vibrating in that mode, that all points on the 
structure will not reach their minimum or maximum values 
Simultaneously. It also indicates that nodal points on the 
Structure will not be stationary. They will have to be 
thought of as instantaneous nodal points. The complication 
of the complex modes can usually be ignored when studying 
lightly damped structures as modal vectors will only stray 


slightly (a 5 or 10 degree phase shift) from being totally 


real valued. 
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If this normalization process is done individually for 
every mode, the relative amplitudes of the modes will not 
have been taken into account. Compared to the lowest 
frequency modes, higher order modes typically have much 
smaller dynamic deflections, due to larger internal damping 
of the more complicated mode shape. Calculating the relative 
deflections of the modes in a Structure can be done by 
comparing the modal parameters between modes on a point by 
point basis. If the 'inertance form'(acceleration/force) of 
the transfer function was used in the modal parameter 
extraction, it must be remembered to divide the residue 
values by a factor of '(2maf)*' to convert these numbers to 
their equivalent ‘dynamic compliance' representation, from 


which the amplitudes can be derived. 


_ re eer 


OTN iin ead 3 ag 2! geescsg Ho ab: 
“$08 TLt¥ Bebo, ec4 3c aipadichit 
: +eeuct nett ot Beregine} FOS 
noavm eves qilesiey: s2000 ebaliae 
patcush fecrerns TOSTH oF: mas ang ba 9 Diab ue 
svitelss an¢d piideivealeo .segare adem berenblgyes 
¥o. @H05- ad aes St s ad eghoer. ads to 
¥d 7015q.'6 Ao 2aGom Ase«l sq e792 9h6tNq Cakigue sa pata 
to (aos0) \46) 347 229008) lh 3 EANAs font! acts o£ : aiesd a 
jesemetean Leon eric « - asf ew fed Samar’ i 
Subse SED) Shiv) of hi 3 acl Ttavn t¢ 4h 


Sow of 


4. PRACTICAL CONSIDERATIONS 
Implementation of the transfer function method of testing, 
involves the simultaneous recording and analyzing of the 
forcing and response signals between designated pairs of 
points on the structure. This chapter deals with the three 
aspects of testing that were considered before and during 
the actual test measurements. The first section deals with 
the selection of the type of loading that was used to force 
the structures. This is then followed by a discussion of how 
field data was recorded employing commonly available 
scientific equipment. The final section deals with coherence 


and its role in verifying the test measurements. 


4.1 Choice of Loading 
There are many types of loading that may be applied 
when making transfer function measurements. These methods 
fall into three major categories. 
1) Sine wave excitation 
2) Random loading 
3) Transient testing 
The use of Sine wave testing has been used successfully 
for over 30 years. Among its early uses were applications in 
the aerospace industry for determining modal parameters of 
aircraft wings. Sinusoidal testing is typically applied in 
one of two methods. In the first, one frequency is excited 
at a time and the response is measured once steady-state 


conditions have been reached to determine the gain and phase 
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shift. The alternative method is a swept sine approach, 
where the structure is forced through a certain frequency 
range and the response to this input is simultaneously 
recorded. Sinusoidal testing offers the advantages of only 
requiring analog equipment and the ability of taking 
accurate meaSurements because the excitation energy can be 
confined to as small a frequency band as desired. Its 
disadvantages include; long testing times, problems with low 
frequency excitation, and the need to set-up, align and use 
Shaker tables (costly and cumbersome). 

Random loading is capable of broad frequency range 
excitation using random signals to drive shaker tables. If 
desired, this random signal may be band-pass filtered to 
allow controlled excitation as was the case with swept sine 
loading. Although random testing Similarly suffers from the 
need to use shaker tables and the fact that windowing must 
be used because it is not periodic in the time record’, it 
does offer the advantage that averaging can be used to 
reduce non-linear effects in the structure and noise in the 
Signal. 

Transient testing, and in particular impact loading, 
offer a simple method to determine the response of a 
structure. Loading is commonly accomplished with an 
instrumented hammer (a hammer having a load cell mounted 
inside the head) which makes it very portable and easy to 
employ for field applications. The impulsive nature of the 


>some signal generators produce 'periodic-random' signals 
that alleviate this problem 
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forcing function produces a broadband excitation signal and 
is to some extent controllable by changing the 
characteristics of the hammer head. Problems with this type 
of loading include: low energy densities in the upper 
frequency ranges, problems with non-linear response of the 
Structure and noisy measurements due to the low signal to 
noise ratio as the signal dies out. Despite these drawbacks, 
the impulse testing method was chosen because of its 
Simplification of the testing process, both in the amount of 
equipment required and amount of time needed to set up and 


uses it 


4.2 Selection of Transducers and Calibration 

Before any field tests were done, it was necessary to 
consider several factors. The size and sensitivity of the 
hammer and response accelerometer were chosen in light of 
the nature of the problem to be studied. For example, a 
light hammer (200 gm.) and a medium sensitivity (10 mV/g) 
pick-up was chosen for the beam test, while a massive hammer 
(5.5 kg.) and a sensitive accelerometer (80 mV/g) was 
required for the compressor base test (low frequency). 

In the tests that were conducted, it was only desired 
to determine the natural frequencies, modal damping ratios 
and relative mode shapes of the structure in question. It 
was therefore not necessary to calibrate the hammer and 
accelerometer. Generally however, this calibration would be 


done before tests were carried out to allow one the option 
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of performing a detailed quantitative study of the test 
results if desired (e.g. actual dynamic deflections or 
determining the mass, damping and stiffness matrices from 


the derived modal parameters). 


4.3 Test Point Determination 

The complexity and size of the structure dictates the 
number of points on the structure to be tested. For the 
free-free beam, thirteen test points were chosen on one-half 
of the beam to get a good definition of the mode shapes. 
Here, use was made of the symmetric nature of the problem to 
reduce the number of test locations required. For the 
compressor base test, again only a minimum number of test 
points were necessary because of the rigid body nature of 
the vibration. Extra test points were included to test the 
validity of the mode shape assumptions. These measurement 
points were chosen so as not to include any suspected nodal 
points of the modes in the frequency range of interest. In 
this way, the derived transfer functions all contained the 
Same number of modes and better global averages of the 
damping and natural frequency could be made. The test points 
were well marked to facilitate a consistent impact location 


and their position on the structure was recorded. 
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4.4 Equipment and Testing 

Figure 4.1 shows a schematic diagram of the set-up for 
recording test data on tape. 

The outputs of the instrumented hammer and response 
accelerometer were first fed into their own signal 
conditioners and then into an amplifier set to boost the 
Signals to fill up as much of the usable voltage range of 
the tape recorder as possible (1 volt RMS). To further 
increase the Signal to noise ratio, the signals could be 
passed through low-pass filters (not shown). In this way, 
only signals in the frequency range of interest would be 
recorded. This would increase the Signal to noise ratio ina 
given frequency range, however it would not allow further 
analysis at a higher frequency range without taking more 
measurements. Most FM tape recorders have a Signal to noise 
ratio of about 40 dB which limits the dynamic range of the 
measurements. As well, very high taping speeds must be used 
with these recorders if it is desired to capture high 
frequency signals (in the recorder used, a tape speed of 60 
in/sec was necessary for a flat frequency response to 207000 
Hz). If low frequency information is not of interest (i.e. 
that below 20 or 30 Hz) direct recording can be used. It 
offers the advantages of higher signal to noise ratios and 
excellent frequency response at low tape speeds (50 to 
19,000 Hz @ 3 3/4 in/sec). 

To ensure that only good signals were recorded, a two 


channel storage oscilloscope was used to monitor the tests. 


to}: 


i 


a ; 4 . i - ia 
OOP. 0S oe stachess qaasugess SEIS Ga cved il etal acw seen - 


a. i-gestezat to don, 6) néeeemioss: | oat ale 1 wot @h* rent | 


bis totsa: seien oc Lengie wedpid to ears oneview siinoasito, 


saneges2 Snes ssmmeil Begnsmstent “ay 19 SJugeve or 


. ont ie 


su yoe Ste fo matpei bo. 7 saresine a 3 a. sawp § 
qed 16 bee pat 


“ - 


iOie WO, Biegl Of02 FSc tas + ‘@vew sasamo ts 2 
os 16 Sead hand ite ares lth® 
sidan i>. 29 fou we By bie 62 oteng 


(2vn uc’ ii. sigiggeg a2 seeTenes 


if 

* 
tr 

’ 
? 
y 
_ 
“=< 
° 
id 
{ 


ieapra wis sa 
aie (yiWore wrli? eescqewal Apeosls 
Siucw! i¢@ate tr sOGnes (Sheepe ys sdtont elengre 
aT [> iscbip-ag? segerad, iteew aden Vosiae 

piic Ton >4 x SSTAVOT . APRS2 y saape 2! 
af L ( epost “udeyors i oapia & 36 ereyte 


jy. op Letpha.¢ avert G1da vere" age) fe tebe <asaemetErm 


mo simaty® Cia ode] Modae Gh te seeds to of3 


mel - thevse tnront’ Apes Yiev ,bLev. ch eee 
| . Pre 
\ 7 Ze ; 
ints eregged os -batcseG.ei sii PNAreoet saads itiv 


belees geo) & chee aeiaee etz of) afeopie qans 


= 


ma)! 
¢i-phebu 4d noo parbyoses Suerte 'sH [2 so 0S voind tet, 


627 96) aheegs aged wel-7s sttioges: «sis pEey soellesxe 


isan .nt £Vé C'S anal _ 


ove «2 .€4812c85 et4w Shanbie note Oi oo. 2nd? ae 


ya 


eae gies ‘G:tnon of Bese aay syoured ( \1Dz0 sensei ais 


a | ay 


pe f + anne 


Accelerometer 


FM 
Tape 
Recorder 
Test 
Structure 
Signal Operational 
Conditioner Amplifier Dual Channel 
Storage 
Oscilloscope 
Instrumented 
Hammer 
Accelorometer - PCB Piezotronics Model 302A 
Hammer - PCB Piezotronics Model 208 AO03 
Signal Cond. - PCB Piezotronics Model 480A 
FM Tape Recorder - Ampex Model FR 1300A 
Scope - Hewlett-Packard Mode! 141B 
Operational Amp - Custom 
Figure 4.1 Equipment used to Record Data 
ia Ge) 
F Floppy 
Operational Low Pass A/D Converter Disk 
FM Tape Recorder Amplifier Filter and Computer Storage 
Krohn-Hite DEC MINC-23 
Ampex (Custom) Model 3342 Filter Data Acquisition 
FR 1300A (96 dB/Octave) System 


Figure 4.2 Equipment used to Digitize Signals 
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This is particularly important when tape recording impulse 
Signals. When using a hand-held hammer to excite a 
Structure, It is difficult to maintain a constant force’ from 
one run to the next; therefore some strikes will be too weak 
and have a poor: signal to noise ratio, while others could be 
too large and saturate the tape. Another problem often 
encountered, was that of multiple impacts between the hammer 
and the test structure. This occurred when a mode with a 
large amplitude was excited and the structure rebounded 
against the impactor before it could be removed. This type 
of loading (quasi-periodic) often causes zeroes in the force 
spectrum (i.e. low signal to noise ratios at those 
locations) and therefore the runs were rejected. 

In actual testing a minimum of averaging (4 averages 
per point) was done on each point to take further advantage 
of the time saving that is inherent with this impulse 
loading method. During both the free-free beam and 
compressor base tests, it was found to be important to allow 
some time between each impact to ensure that the structure 
had completely stopped vibrating from the previous 
excitation. If this was not done, the coherence between the 
runs would have been affected as the response measured from 
the ensuing impact would not have been solely due to that 


impact (see discussion on coherence later in this chapter). 
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4.5 Digitization of Signals 

Once field neasurementSiiwere (complete, the irecorded 
Signals were digitized for processing. Figure 4.2 shows the 
arrangement of the equipment used in digitizing the data for 
this study. The output of the tape recorder was first 
amplified to fill the largest amount of the full scale range 
of the analog to digital(A/D) converter as possible (the 
digitizer used had only one voltage range; +5 volts 
maximum). The signals were then passed through anti-aliasing 
filters (sharp roll-off, low pass filters @ 96 dB/octave) 
and then into A/D converters. This data was labelled and put 
on floppy disks for further processing (curve fitting, 
etc.). With the particular digitizer used, the digitization 
was not een Simultaneous between channels. This type of 
Situation could have distorted the original relationship 
between the signals, however owing to the sampling rates 
that were used, this problem was safely ignored (see 


Appendix B). 


4.6 Coherence 

The development of the transfer function relations in 
Chapter 2 assumed that the measured applied forces to the 
structure were its only inputs and that the response was due 
solely to SeHG Situ In practical testing however, the 
recording instruments have their limitations and test 
objects can not be totally isolated from external influences 


so that this condition is not completely met. Typical causes 
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for this include, electrical noise, limited dynamic range of 
the recording instruments, non-linearities in the 
transducers and the structure as well as vibration induced 
from the surroundings (e.g. rotating machinery, busy roads). 
In order to measure the effects of these external 
2 


influences, the coherence function, y*, was used. This 


function can be calculated from the following expression 


WCE) = Gree aee sh) 


Ga.\ fonGey bh) 


where G,,(f), Gy, (f) and G,,(f) are the averaged input and 
Output power spectrums and cross spectrums respectively (the 
Star denotes complex conjugate). The coherence function 
which can only take on values between '0' and '1', can be 
thought of as a measure of the power in the response signal 
caused by the power in the recorded input. If the coherence 
equals one at some frequency, it implies that the output 
power at that frequency was entirely due to the measured 
forcing signal. When the coherence function takes on a value 
other than one, it shows that the system response at these 
frequencies was not totally due to the input. 

Ideally then, the coherence function should equal one 
over the entire frequency range of interest. By graphically 
displaying this function along with the determined transfer 
function, serves aS a valuable aid in determining certain 
properties about the recorded signals. In the tests that 
were done during this study, drops in coherence could 


invariably be attributed to the small dynamic range of the 
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PM tape recorder. Figure @.3 shows a typical transfer 
function obtained from compressor base tests. In this 
diagram, it is seen that wherever the frequency content 
falls approximately 35dB below the peaks in the spectrum, 
that both the transfer and the coherence function are 
affected. In these regions the transfer function is no 
longer a smooth curve and is almost random in nature. This 
is due to the fact that an FM tape recording system has a 
dynamic range of about 35 to 40dB. When signals are recorded 
with a larger dynamic range than this, the smaller amplitude 
Signals will not be recorded. This also causes the coherence 
to become very jagged, because in these frequency ranges the 
true signals were not recorded and what remains is a noise 
floor from the tape. 

In the tests conducted, it was not necessary to have a 
large dynamic range so that the recording means was 
Satisfactory. If a larger dynamic range is required several 
options can be considered. First, the tape recording aspect 
of the field data can be totally eliminated by directly 
digitizing the signals on location. Employing a twelve-bit 
A/D converter, a 72dB dynamic range can be realized. 
Secondly, if the flexibility of tape recording is desired 
and low frequency signals are not expected in the results, 
direct recording can be used. Typical direct (AC) recording 
systems have dynamic ranges in excess of 70dB and have large 
frequency bandwidths. When low frequency tape recording is 


necessary, a thirdwaiternative: is the usevof an FM tape 
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recording system that employs a switching unit to record one 
Signal on two channels of the tape. One channel records that 
portion of the signal that lies in the zero to -35dB range 
while the second channel records anything below the -35dB 


level. 
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Figure 4.3 Coherence of Tape Recorded Signal 
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5. TEST RESULTS 

In the two previous chapters, the procedure for 
determining the modal parameters as well as an outline for 
collecting experimental field data was presented. In this 
chapter the details of tests conducted on a free-free beam 
to verify the developed procedure are given. As well, 
results of an application of this technique to the dynamic 
analysis of an air compressor and its inertia base are 


included. 


5.1 Free-Free Beam 

A free-free beam was chosen for testing because of the 
potential ease of simulating the boundary conditions. A 
cantilever beam was considered, but in preliminary testing 
it was found that the measured response and corresponding 
frequency spectrum contained extraneous signals that could 
be attributed to the 'built in' support. This would not have 
caused any problems in the analysis of the derived transfer 
functions, but it was desired to do the first modal test on 
a structure that gave the cleanest Signal possible. 

A hot-rolled steel bar of dimensions 9.46 mm. x 37.4 
mm. x 539.3 mm. was used as the test beam. Calculations 
showed that three natural frequencies for this bar 
(excluding rigid body modes) existed below 1000 Hz. Three 
widely spaced natural frequencies with low damping, was 
deemed a good starting test for the analysis system. The 


free end conditions for the bar were simulated in early 
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tests by placing the bar on a foam support. Two major 
problems were encountered using this technique. The first 
problem was that it was very difficult to impact the beam 
and pull the hammer away before the beam would bounce back 
and again make contact with the hammer. The discussion in 
the previous chapter pointed out why this type of signal was 
a problem. The second problem encountered using the foam 
base, was the very high damping induced by the support. In 
the analysis, this effect showed up as high modal damping 
parameters which were not expected in a steel structure. To 
overcome both of these problems, the bar was suspended by 
elastic bands at both ends. These low frequency supports 
induced rigid body modes that had natural frequencies in the 
order of 2 to 3 Hz; well below the frequency range of 
interest. 

Because the study was to be restricted to those modes 
below 1000 Hz, it was not desirable to put any energy into 
the higher frequency ranges that were not to be analyzed, 
Since this would lower the signal to noise ratio in the 
testing range. To meet this requirement, a hard rubber 
impact tip was used with the hammer. The resulting input 
signals had a spectrum that was flat (+3dB) to approximately 
1300 Hz and decayed gradually beyond this range. 

From Chapter two, it was seen that it is necessary to 
define one row or column of the transfer matrix to be able 
to completely describe the dynamic behavior of a structure. 


This is equivalent to taking measurements with a fixed 
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response pick-up and moving the force input to the desired 
locations on the structure (row determination) or fixing the 
forcing location and changing the position of the output 
transducer. (column.determination). In this test,.it.was 
chosen to rigidly mount the accelerometer at one end of the 
beam and impact the structure in different locations along 
its length because of the ease in moving the hammer from 
position to position. Only one half of the beam was tested 
due to the symmetry of the bar and itsS support conditions. 
By mounting the accelerometer at the end of the beam, it was 
assured that the response of ae beam would contain 
components of the lowest three modes in question, as all 
three have large modal deflections in this region. If the 
accelerometer had been placed at a node of one of the modes 
that it was desired to analyze, the response of that mode 
would not have been picked up, irrespective of where the 
beam was excited. 

To choose the forcing locations along the beam, a 
single channel FFT analyzer was used to monitor the response 
to different impact locations. At locations where all three 
modes were present, a piece of tape was used to mark the 
location. If a nodal point for one of the modes was located, 
positions were tried on either side of the Original spot 
untadls as suitable’ areaswas.found.-In thesend; antotal.of 
thirteen impact locations were located on the one half of 
the beam at approximately 20 millimeter intervals. (It is 


important to note here that, this trial and error type of 
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approach for finding suitable test locations is somewhat 
Slow and cumbersome and resulted from the lack of foresight 
in the original program development that analyzed the 
experimentally determined transfer functions. Generally 
however, testing points could be chosen arbitrarily and the 
analysis programs employed would account for the fact that a 
particular mode of vibration might have a nodal point at the 
location in question.) 

AS waS mentioned in the previous chapter, an FM tape 
recorder was used to record the input and response signals. 
In these particular tests however, a direct recording unit 
would have been sufficient as there was no low frequency 
information (below 20 Hz) that was of interest. 

Table 5.1 lists the averaged damping ratios and 
resonant frequencies determined from the tests and compares 
these resonant frequencies to those predicted theoretically 
for a free-free beam. The low damping ratios correspond with 
those typically observed in steel structures. For all three 
modes, the predicted frequencies are larger than the 
experimentally determined numbers. Comparing the ratios of 
f,/f, however, very good agreement is seen. This indicates 
that the two cases differ by a scale factor. This scale 
factor may simply be due to the difference in the material 
constants (Young's modulus and mass density) and geometry of 
the bar compared with those assumed in the theoretical 


calculations. 
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Table 5.1 Experimental vs. Theoretical Resonant 
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The experimentally derived, normalized mode shapes are 
given in table 5.2. Examining the phase angles of the modes, 
it is seen that parameters very nearly describe normal 
modes; that is, the modal vectors are either in phase with 
each other (phase angle near zero), or directly out of phase 
with one another (absolute value of the phase angle close to 
180°). Figure 5.2, compares the experimentally determined 
mode shape (disregarding the imaginary component of the 
normalized modal vectors) to that given from the solution of 
the analytical equation. It is seen that the experimentally 
determined mode shapes agree very closely with their 
predicted behavior. (It should be noted that all three mode 
shapes have been normalized relative to the corresponding 
modal parameters of the first point on the beam. In doing 
so, the relative amplitudes of the modes have not been 


considered). 


5.2 Compressor and Base 

The second set of tests were performed on an 'Atlas 
Copco' air compressor and inertia base system that was 
designed to run at 720 rpm(12 Hz). The air compressor was 
driven by an electric motor with a pulley drive. The 
motor-pulley system and air compressor were mounted on a 
steel base filled with iron shot. The total mass of this 
system was approximately 3100 kg. with 580 kg. due to the 
compressor alone. The configuration was isolated from the 


floor by 4 steel springs, each with a stiffness of 
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Figure 5.1 Mode Shape of the 167 Hz Mode 
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approximately 670 kN/m. 

The simplest model of this setup would be that of a 
rigid body, free to move in all directions. This model would 
imply that six rigid body modes should be present in the 
structure. If the mass and stiffnesses were symmetric in all 
three planes, these modes would be uncoupled as shown in 
migure (5.4. 

Any non-symmetries in the structure would tend to 
couple the modes in such a way that a mode would no longer 
be purely translational or rotational, but rather some 
combination of the three translational and three rotational 
motions. Because of the potential coupling, it is possible 
that all six modes could have been excited from one loading 
direction. However, small amplitudes resulting from any 
weakly coupled modes could have made the identification of 
reliable modal parameters difficult to obtain. It is also 
quite possible that some modes would not have been excited 
at all in the particular direction of loading considered and 
therefore not detected in the results. With this in mind, it 
was chosen to excite the body from three different test 
locations; one in each of the three principal planes. In 
this way the chances were good that all of the modes would 
be excited from at least one of the forcing locations 
(unless the three loading points happened to lie on a node 
of one the modes in question). 

The testing points on the compressor system are shown 


in’ figure: 5.5 ..Thesthiokwarrows (labelled A™, '5' ‘and \C') 
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indicate the impact locations while the other arrows show 
the points where the response was measured. The test 
locations were chosen as far apart as possible to give a 
better definition of the mode shapes and yet in a symmetric 
manner to simplify the analysis. 

In the study of the free-free beam, the accelerometer 
had been held fixed and the hammer moved from point to 
point. Due to the layout of the compressor base system 
however, it was chosen to have a fixed forcing point in each 
testing plane. This was done because the large hammer used, 
made it difficult to ensure good clean impacts unless the 
location was easily accessible and had a flat solid surface 
(i.e. it was desired to hit the surface squarely and not 
damage any of the structure in the testing). 

Preliminary tests on the structure determined that all 
of the modes of interest were below 10 Hz. The hammer that 
was to be used however had a much broader excitation range; 
the softest head available for the hammer excited a flat 
frequency response ‘to about 100 Hz. When this signal was 
recorded on tape and later low pass filtered below 10 Hz for 
digitizing, very little of the signal was left to work with. 
To avoid the low signal to noise ratio associated with this, 
it was decided first to further amplify the impulse signal, 
then low-pass filter it before recording to maintain as high 
a signal to noise ratio as possible. In implementing this 
procedure, some of the inherent advantages and flexibility 


of this type of recording system were lost, (i.e. it would 
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be necessary to do further testing if higher frequency 
analysis ranges were desired) but it was felt that this 
trade-off was necessary. 

The following three figures show examples of transfer 
functions obtained from testing in each of the three planes. 
Four distinct modes are present in Figure 5.6 which were 
obtained by loading the base at point 'A' (in figure 5.5) 
and recording the response right beside it (point #1). The 
rounded broad nature of the peaks would indicate that heavy 
damping was present in the system. 

The transfer function formed from forcing the base at 
point 'B' and picking up the response at point #5 is shown 
in figure 5.7. In this diagram, five major peaks are seen 
between the frequency range of 3 Hz and 6.5 Hz. (It was in 
this range that all six rigid body modes were found.) It 
should be noted that three peaks located at 3 Hz, 4.7 Hz and 
5.4 Hz are common between this figure and figure 5.6. 

Shown in figure 5.8 is the transfer function derived 
from forcing the structure at point 'C' and recording the 
response at point #10. Two modes are located near 3.1 Hz (as 
was the case in Figure 5.7) as well as two other modes 
around 6.2 Hz. The presence of two distinct higher frequency 
modes is not obvious in the magnitude plot of the transfer 
function; the phase angle of the function must be studied to 
see that indeed two separate modes exist in this region. 

The modal parameters where then obtained for the six 


modes by curve fitting the measured transfer functions. The 


, 
- . 


ams = : 
ynaspat% ragga a unitess slat Pry 
"ve 
Sees) 7 | hie al ii ‘ae ‘ss 
oe) \ “4 vé 
3 7 see h i . 
ayi te 2stgageehwete casual? >eeney ee 2. 
5 He i = rT 


equelad esd? Sit $6 eee © (teak Oe: Beep 260. anc ay 
sv doinaw €@.@ Sets fy -7ekste << oi BOM 20 se tb. 

vc. 

PQUPiS NEV A* “7 Fe-Oq, a = 8 o2at of9 peaieel qa £ fabbie : 

. 3h2eE ; wb free tity ior Serer, es oni! 7 

a2 4325 *! ow 2 Nee jefe 10! sindiais beard. | 

nS ogee at. godesg. 2a ie 

a2/a5 ! Aanivio: £f 7%, Cea mMmitgnet vsherass 

, 'h'So oo 42neees ‘ota au gniaute Sas ‘é* 


:2 on ; >. nl » ts Sr ey hae Sabet bt Sqn cty At a o1pé . 


if 


cH 2.9 Git oF f Be song ztnwupes’ eas. 

eine bree atom (hed Sapo ate Die aus mia 
zt wites souks tae boson od, 
se. =" Saure 24 a sures git: anew dad homey e>6 3 
| wtih, $atarw: 5, OAS BER ap epbgth. at owert@ 
oti pashstos7 08 = PGR apis Sevdputee | tigi painiad: 
26) 361 2620) Deseret. (=e webioa 97.0% syikog te) | 
Sabha 26A4s Sud oe’ rr 7 Gre espe ot od 
oinaibests =4cn14 Pe ee sonae:4 sd? a $08 
| Helknats of2 35 soig Sess ERR 99 ne eve igen aEts 
eo Reihuca: ad teum italia Sips “en 
ngiees ists: nt sean aatibe: ‘Sraueqre cws 
ute ae. joy bunt sate waa ayecasioes 


. 
* 
in 
a 

oF 
om 
4, 

Lg 


t 


t 
F 
iY 
* 
‘ 


60 


a 9 VAT 
~180.0 pee MUA Gili J We 


FREQ 10.0 
TP # 1.0 


Figure 5.6 Transfer Function in X-Y Plane 
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tightly spaced modes near 3 Hz and also those around 6.2 Hz 
caused inconsistencies when curve fitting because of the 
poor resolution and high overlap of the two modes in either 
region. This was further hampered by the exponential window 
(which increases the apparent damping) that was used on the 
response data to help reduce the noise in the results. (It 
should be noted that the damping in the system was very 
heavy; the response signal decayed completely in about 
one-quarter of the time window, leaving the majority of the 
Signal with a high noise content.) 

Table 5.3 lists the averaged resonant frequencies and 
damping ratios that were obtained from the curve fit of the 
rigid body modes. The large damping ratios obtained, 
illustrate the heavy damping that was present in the system. 

Shown in figure 5.9 is the mode shape that occurred at 
a resonant frequency of 6.51 Hz (only the base of the system 
is shown to simplify the representation of the mode). This 
mode shape displayed an almost purely rotational mode about 
the x-axis. The mode shape that was associated with the 6.30 
Hz mode is shown in figure 5.10. This was mainly a 
rotational mode about the z-axis with some light coupling 
from other modes (especially the rotational mode about the 
Y-axis). 

The mode shape associated with the 5.44 Hz resonance 
was unreliable because of problems associated with curve 
fitting. What was clear from the results however, was that 


its major component was that due to rotational motion about 
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MODE DAMPING RATIO RESONANT FREQUENCY 

(%) (Hz) 
1 T53 S06 
2 0.80 See 
s C25 4oi2 
4 Le 30 5.44 
5 1750 6.130 
6 O55 6.51 


Table 5.3 Resonant Frequencies and Damping 
Ratios of the Compressor System 


a line parallel to the y-axis. This behavior is illustrated 
in Figure 5.11. More coupling from other modes was probably 
present than what is depicted in this diagram because this 
rotational motion was also seen in the 6.30 Hz, 4.72 Hz and 
3.06 Hz modes. The reasoning behind this is that coupling 
with the other modes would have been necessary if the effect 
of this particular motion was viewed in other modes as well 
(i.e. if two modes are coupled, the effect is seen in both 
resulting mode shapes). 

The mode shape of the 4.72 Hz peak is shown in figure 
5.12. This mode consists mainly of a vertical translation 
with some coupling from the 5.44 Hz mode to produce the 
Small rotation about the y-axis. 

Thus far, the modes the have been defined have included 
three mode shapes that were primarily rotational in nature 
as well as a fourth mode that had a mainly vertical 


translation motion. Of the six expected rigid body modes, 
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this leaves only the two horizontal translation modes to be 
identified. The mode shapes determined for the resonant 
Brequencres of 3.06 Hasand3,22 Hz are shown if figures 5713 
and 5.14 respectively (note - the close proximity and heavy 
damping of the two modes produced inconsistent modal 
parameters, so that the mode shapes shown are estimates 
only). Both mode shapes primarily displayed a translational 
motion with the 3.06 Hz mode translating mainly in the '‘y' 
direction, while the 3.22 Hz mode had the largest component 
of its displacement in the 'x' direction. If the mass 
distribution and lateral spring stiffness of the structure 
had been perfectly symmetric, it could be shown that two 
purely translational modes with identical natural 
frequencies would exist. These modes would exhibit 
physically perpendicular mode shapes with no preferred 
orientation of movement. When either the spring stiffnesses 
or the center of mass location is non-symmetric however, the 
two frequencies become distinguishable from one another, the 
modes have a fixed orientation and the mode shapes no longer 
exhibit perpendicular motion. Therefore in examining the 
results for the two translational modes, it is not 
Sunprising to find that thers resonant frequencies are so 
close together or that their mode shapes show some motion in 


both the x “and Myo rdirections. 
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Figure 5.9 Mode Shape of the 6.51 Hz Mode 


Figure 5.10 Mode Shape of the 6.30 Hz Mode 
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Pigure 5.13 


Mode Shape of the 5.44 Hz Mode 


Figure 5.12 


Mode Shape of the 4.72 Hz Mode 
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Figure 5.13 


Figure 5.14 


Mode Shape of the 3.22 Hz Mode 


Mode Shape of the 3.06 Hz Mode 
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6. CONCLUSIONS 

The testing and analysis regime that has been presented 
in this study has focused on the practical implementation of 
the modal analysis technique. Developed during the course of 
this study was an economical data acquisition system and 
analysis software that could be readily adapted by smaller 
engineering firms that could not justify the expense of a 
sporadically used all-in-one instrumentation system. 

From the work done in the study, two things are clear: 
the first is that good reliable transfer function data can 
be recorded using commonly available instrumentation. (In 
fact, the use of a portable tape recorder may be deemed 
essential in testing locations where on site digitization of 
the signals is impractical). The second thing to be noted is 
that modal analysis is not a 'blind' technique that may be 
effectively applied with little prior knowledge about the 
method. Familiarity with this type of testing is essential 
in taking good field measurements and later properly 
interpreting the derived modal parameters. 

The modal parameters identified from the transfer 
functions of the free-free beam were consistent and 
reliable. However in dealing with the compressor system, 
many problems associated with parameter identification were 
encountered. For the most part, the difficulties could be 
attributed to the close spacing of some of the modes and the 
heavy damping in the system. Despite this, it was still 


possible to separate and identify the modes in question. 
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However if greater accuracy and more consistent results were 
desired, it would be necessary to first of all excite the 
Structure in such a way that more energy was being 
concentrated in the lower frequency range and also to employ 
‘zoom' Fourier transform techniques to allow greater 
frequency resolution in the range of interest. 

In the course of the study, the area that required the 
most research time was that involved in curve fitting. The 
iterative Levy approach that was used in the initial 
analysis of the transfer functions proved to be very 
reliable in obtaining preliminary parameter estimates. When 
further curve fitting was necessary to reduce the errors 
between the experimental and analytical functions, a good, 
quick method was not found. Further study in this area would 
be necessary to research and develop software appropriate 
for smaller computing facilities. 

The next logical stage of development of this modal 
analysis testing system would be to make use of the modal 
parameters derived from curve fitting to define a lumped 
parameter model for the structure being studied (by 
calculating the mass, damping and stiffness matrices of the 
system). This would then enable the possibility of using 
structural modification techniques to quickly determine the 
effects of any proposed changes to the structure. This type 
of analysis would be useful in the retro-fitting of existing 
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APPENDIX A - LEVY PROGRAM 


SUBROUTINE LEVY (HEXP,FREQ,NLOW,NHIGH,K, VALUE) 


» COMPUTES A LEAST-SQUARES FIT FOR A TRANSFER FUNCTION 
OF UP TO FIVE MODES USING LEVY'S METHOD 


HEXP,FREQ - EXPERIMENTAL COMPLEX TRANSFER FUNCTION 


AND CORRESPONDING FREQUENCY LOCATION 


NLOW,NHIGH - LOWER AND UPPER LIMIT OF CURVE FIT RANGE 
Kia, NUMBER NOP POINTS TO BE FIT IN THE RANGE 


NOTE 


REAL*8 
REAL*8 
REAL*8 
REAL*8 
REAL#8 


-- VALUE - OUTPUT MODAL PARAMETER VECTOR 


CALCULATIONS WERE DONE IN DOUBLE PRECISION 
BECAUSE 'C' MATRIX WAS FREQUENTLY 
POORLY CONDITIONED 


C(20,20)/400*0.0D0/,D(20)/20#0.0D0/,WKAREA (500) 
RE(401),IM(401),SUM(401) ,VALUE(20),TWOPI ,FF 
ABSDEN(401)/401*1.0D0/,TEST,B1,B2,F(21,401) 
C10) 710*0.0D0/7 -TO10) 7 10*0..0D0/,S( 10) 710*0..0D0/ 
U(10)/10*0.0D0/,A(11),FREQ#4(401) 


COMPLEX#16 Z(10),DEN,ALPHA(10),AA,BB,HEXP#*8(401) 
COMPLEX*16 ZOLD(10)/10*(0.0D0,0.0D0)/ 
TWOPI=2.0D0#3.141592653589793 


K2=K*2 
K4=K*4 


K4P1=K4+1 


PUT POWERS OF FREQUENCY INTO MATRIX F(I,J) 


DO 100 


F(2,J)= 


J=NLOW, NHIGH 
FREQ(J) 


RE(J)=REAL(HEXP(J) ) 
IM(J)=AIMAG(HEXP(J) ) 
TREDABSGR(2 3 )9ET..0.0000001) F(2 J) =0 70000001 


F(2,J)= 


F(2,J)*TWOPI 


IF(K.GE.3) F(2,J)=F(2,J3)/1000.0D0 


SUM(J)=RE(J)*RE(J)+IM(J)*IM(J) 
DO 100 I=1,K4P1 

Fir Wi ]P (270) 52 Gr) } 

CONTINUE 

IC=0 

GOfTo+104 


CORRECTING FACTORS FOR REDUCING WEIGHTING EFFECTS 


DO 103 J=NLOW,NHIGH 
AA=DCMPLX(0.0D0,F(2,J)) 

DEN=DCMPLX(1.0D0,0.0D0) 

DOMMO Zw 1.7K 

DEN=DEN+D(K2+1*2-1) *AA¥*(1*¥2-1)+D(K2+1*2) *AA¥* (1 *2) 
CONTINUE 
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ABSDEN (J) =DEN*DCONJG (DEN) 
103 CONTINUE 


COMPUTE CONSTANTS THAT GO INTO 'C' AND 'D' (CX=D) 


ANAND 


104.5DG? 105321 -K2 
GO )=070D0 
Tio) =O 0 D0 
S(J)=0.0D0 
U(J)=0.0D0 
105 CONTINUE 
DO= 106 F=1)7K2 
DO 106 J=NLOW,NHIGH 
L(I)=L(1)+F(2*1I-1,7) /ABSDEN(J) 
T(I)=T(1)+F(2#I ,J)*IM(J) /ABSDEN (J) 
S(1)=S(1)+F(2*I-1,J)*RE(J) /ABSDEN(J) 
U(I)=U(1)+F(2*1+1,7)*SUM(J) /ABSDEN (J) 
106 CONTINUE 


ELEMENTS IN THE VECTOR 'D' 


GOO 


DO 107 I=1,K 

D(2*I-1)=S(I) 

D2) =T C1} 

D(2#*#(K+1I)-1)=0.0D0 

D(2*(K +I ))=0 bE) 
107 CONTINUE 


ELEMENTS IN THE MATRIX 'C' 


AAA 


DOM1GS “= b/K 
NEG=(-1) **(J+1) 
Poelgo. 1= 1K 
I2=1*2 
121=I12-1 
G2 ae 2 
J21=J2-1 
TJ=T+J 
IJ1=IJ-1 
C(121,021)=L(IJ1)#NEG 
C(K2+121,J321)=T(IJ1)*#NEG 
C(K2+12,J321)=S(1J) *NEG 
C(12,J32)=L(1J) #NEG 
C(K2FI2T 32) =S (1d )'* (-NEG) 
C(K2+12,J32)=T(1IJ) #*NEG 
C(I21,K2+J21)=T(1J1)*NEG 
C(12,K2+J21)=S(IJ)*(-NEG) 
C(K2+121,K2+J21)=U(IJ1)*NEG 
C(121,K2+J2)=S(1J) *NEG 
C(I2,K2+J2)=T(IJ)*NEG 
C(K2+12,K2+J32)=U(1J) *NEG 

108 CONTINUE 
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J1i=1 
J20=20 
J0=0 


THE 'IMSL' SUBROUTINE 'LEQT2F' SOLVES A SYSTEM OF 
SIMULTANEOUS EQUATIONS 


GALH LEOTIF (GC) Jit) K41920,D 7/50, WKARBA, IER) 
DETERMINE THE ZEROES OF THE DENOMINATOR 


DOMA9! L=ToK2 
A(1I)=D(K4P1-1) 
CONTINUE 

A(K2+1)=1.0D0 


THE 'IMSL' SUBROUTINE 'ZRPOLY' COMPUTES THE ZEROES 
OF A POLYNOMIAL 


CALL ZRPOLY(A,K2,2Z,IER) 
CHECK RELATIVE ERROR OF DAMPING AND NATURAL FREQUENCY 


TEST=0.0D0 
DOehTOT=27.K2 
BB=Z(I)-ZOLD( 
B1=DREAL(Z(1I) 
B2=DIMAG(Z(I)) 

TEST= TEST+DABS (DREAL (BB) /B1)+DABS (DIMAG (BB) /B2) 
ZOLD(I)=7 (Cr 

CONTINUE 

IC=IC+1 

LEXCLC.EO}20)) GOuTO. 111 

IF (TEST/DFLOAT(K2).GT.0.001D0) GO TO 101 

GOnTO 1:13 

WRITE(10,112) 

FORMAT(1X,'LEVY DIDN''T CONVERGE IN 20 ITERATIONS’ ) 


COMPUTE COMPLEX RESIDUES 


DOwiiS ba teR2 
ALPHA(I)=0.0D0 

DO; 114 d= 

Jeade2 
ALPHA(I)=ALPHA(I)+D(J2-1)*2(1)**(J2-2)+ 
#D (02) 42 (1) e*8(a2-19 

CONTINUE 

ALPHA (I )=ALPHA(I ) /D(K4) 

DO - MS Jake 

IF(J.NE.1) ALPHA(1I)=ALPHA(1I)/(2Z(1)-Z(J) ) 
CONTINUE 


CONVERT MODAL PARAMETERS 
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DOr116. l=1,K 

IIl=2*I-1 

VALUE (4*I-3)=DIMAG(Z(II)) 
VALUE ( 4*1-2)=DREAL(Z(II)) 
VALUE (4*I-1)=DREAL(ALPHA(II) ) 
VALUE (4*I )=DIMAG(ALPHA(ITI) ) 
CONTINUE 

PCRs T 63 ts GOe TO. e118 

DOsdided =e Ks 

VALUE(I )=VALUE(I)*1000.0D0 


/TWOPI 


7 CONTINUE 
8 RETURN 


END 


me 


. Son 


APPENDIX B - NON-SIMULTANEOUS SAMPLING 


When using dual or multichannel digitizers, it is often 
assumed that all of the channels are being digitized 
Simultaneously. This is usually not the case however. Most 
A/D converters have some finite time lag between the 
sampling of one channel to the sampling of the second and 
ensuing channels. If this time-lag, t>, is very small 
compared to the sampling period, At, (ie. t o/At > 0) the 
data collected will be valid. As the digitizing speed is 
increased (to is fixed), the relationship between the 
Signals on two separate channels will become noticably 
different from two which are simultaneously sampled. The 
time-shifting property of the Fourier transform pair can be 
used to examine this effect. 

Consider a signal h(t) and its Fourier transform H(f). 
If h(t) is shifted along the 't' axis by an amount to, it 
can be shown that the Fourier transform of this function, 
h(t-t.) is H(f£)exp(-jwt,.). To get the desired frequency 
spectrum, H(f), that would have occurred had the sampling 
been truly simultaneous (to=0), it is necessary to simply 
divide the calculated Fourier transform by the factor 
exp(-jwt,). If the product wt, (which is proportional to 
t,/At) is small, then this factor approaches the value of 
one (real valued only). Thus small sampling delays cause 
little distortion tomthe spectrum of interest. However, if 
the term wt>o, has some appreciable value, then the factor 


exp(-jwt,) will be complex valued and will affect the phase 
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angle of the spectrum (the magnitude stays the same). 

In practical testing then, it is necessary to determine 
what time delay (if any) is present between the channels of 
the digitizer. If the ratio t,./At is small, it may be chosen 
to neglect the correction term (typically if t,. is less than 
1% of At, the correction is ignored). Generally however, the 
addition of the correction factor can be simply implemented 
into the analysis procedure, so that any time delays can be 


corrected. 
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